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The results of an earlier paper, describing pulse propagation in multi- 
mode dielectric waveguides with random coupling, are specialized to the 
two-mode case. Because of their greater simplicity, the results for this 
special case provide more insight into the mechanism of pulse shortening 
due to mode coupling. The two-mode theory yields a formula for the width 
of a pulse carried by coupled guided modes that is found to hold also for 
four modes, so that it may be true for an arbitrary number of modes. This 
formula [eq. (23)] contains only the measurable distance required to 
establish the steady-state power distribution and the length of uncoupled 
pulses. The pulse length formula is identical with Personick's important 
result. Our treatment suggests that the characteristic length appearing in 
this formula may be accessible to measurement. 

I. INTRODUCTION 

In a series of earlier papers, the theory of multimode propagation in 
dielectric waveguides was analyzed with the help of stochastic coupled 
power equations. 1-3 The propagation of Gaussian-shaped pulses in a 
waveguide with randomly coupled modes was treated quite generally 
for N modes in Ref. 3. The theory was based on the following form of 
the stochastic coupled power equations: 

ap i ap N 
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P„ is the average power in mode v, v, the group velocity, a, is the atten- 
uation coefficient of mode v in the absence of coupling to other guided 
modes, and h,„ is the power coupling coefficient. To second order of 
perturbation theory, and assuming a Gaussian shape (in time) of the 
input pulse, the solution of (1) can be expressed as: 3 

P*.0-§^W«V-«-«p{-(^)]. (2) 
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Mi , the width of the z'th Gaussian function in (2), is given by 

A/, = 2(t 2 + faj"*)* (3) 

The input pulse with half width r and amplitude G, is assumed to be 

P, - G, exp (-4)- (4) 

The coefficient fc, appearing in (2) is determined by the input pulse 

k< = tfl#. (5) 

r-l 

The vectors with components J5j° and the parameters «J° are the 
t'th eigenvectors and eigenvalues of an algebraic eigenvalue problem 
denned in Refs. 2 and 3. The parameter a 2 {) is the second-order pertur- 
bation of the eigenvalue «J° and is denned as follows: 



a u-> = ^ U-, ^ g/ . L. (6) 

^t a„ - a 

v is the average group velocity. 

This approximate theory of pulse propagation in multimode wave- 
guides holds for random coupling between the guided modes under the 
assumption that the correlation length of the coupling function is short 
compared to the distance over which the mode power P, changes 
appreciably. 

II. APPLICATION TO THE TWO-MODE CASE 

In its full generality, the theory of multimode pulse operation is hard 
to evaluate. Computer solutions are being provided in Refs. 2 and 3. 
Here we want to derive expressions for the special case of two modes, 
that allows us to gain more insight into the meaning of the theory. The 
two-mode case has been treated previously by several authors. 4 ' 5 Our 
results are thus not all new. 

For two modes, we write h 12 = h 2l = h. We now have to solve the 
eigenvalue problem 

(a„ - a, - h)B l + hB 2 = Ol (7) 

hB x + (a. - a 2 - h)B 2 = 0) 
The equation system (7) has the solution 
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= h + ^±^ _ [fa""') 2 + tfT (8) 



(8) 



ai + a 2 , |~ ( g 2 — «i) 2 , ,2 I* 



(9) 



for the first and second eigenvalue. The two components of the first 
eigenvector are 

B[ " = [53 + ,]•{_ _ i ; 2 [fe^ + ;/]'}' (10) 

ai" = * r L * -,, J J • (id 

2[ fa ~ ai) + A'] 

The components of the second eigenvector can be expressed in terms of 
the components of the first eigenvector. 

B[ 2) =B a " 2? a 2) = -51". (12) 

III. DISCUSSION OF THE TWO-MODE CASE 

In the special case a x = a 2 = we have the eigenvalues 

a[ l) = (13) 

and 

aT = 2h, (14) 

while the eigenvectors are 

B!» = i *»-£ (15) 



and 



B > ] =V2 Bi " = -k- (16) 



There are several interesting features apparent in this special solution. 
In the absence of loss, both modes carry equal power. We see im- 
mediately from (15) and (16) that the sum of the squares of the com- 
ponents of each eigenvector adds up to unity, while the inner product 
of the two vectors vanishes. This is a general property that is also shared 
by the solutions (10) through (12). 
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The "loss coefficient" a"' of the term with i = 1 in (2) vanishes, so 
that this term does not decrease in amplitude as the pulse moves along 
the 2-axis. The "loss coefficient" af of the second term of (2) (with 
i = 2) is equal to twice the coupling coefficient h, so that this term 
becomes vanishingly small for large values of z. This too is a general 
feature of (2). The lowest order eigenvalue a" 5 is smaller than all other 
eigenvalues, so that only the first term of the series (2) remains for large 
values of z while all the other terms have become vanishingly small. 
Even though «"' is not zero in the general (lossy) case, the multimode 
waveguide always reaches a steady state which is described by the first 
term of the series expansion in (2), provided that the modes are coupled. 
Only in the lossless case do we find a^ = 0. It is noteworthy that the 
"loss terms" exp (— «J°*) (i > 1) all become vanishingly small even 
in the absence of losses. A steady-state distribution of mode power 
versus mode number is thus established that is independent of the 
initial excitation of the waveguide. The decay of the higher-order terms 
in (2) does not necessarily indicate power loss. We see indeed, from the 
solution (16), that the sum of the components of the second eigenvector 
adds up to zero indicating that no power is carried by the second term 
(t = 2) of (2) in the absence of loss. The individual terms of (2) must not 
be confused with waveguide modes. They have no independent physical 
meaning except for the first term with i - 1, which is the steady-state 
power distribution. It is apparent that it is sufficient to study the 
behavior of the first term in (2) alone, since all other terms (the second 
term is the only other term in the two-mode case) become negligible 
for large values of z. 

We can easily define a characteristic distance that is required for 
the steady state to establish itself. Once the exponential factor 
exp (a l l) — ai 2) )z has become small, the steady state is reached. We 
thus define the characteristic length as follows 

L > = a < 2 > _ a <» " ' 17 ' 

The parameter k is a number of order unity. For k = 1 we have 
exp (al l) - a 2) )L, = l/e. Thus k = 1 is too small to consider the steady 
state as reached. However, we can still define L. by (17) with k = 1. If 
we use k = 4.6, we have exp (a' 1 ' - a™)L t = 0.01. This number is 
small enough to consider the second term in (2) as negligibly small. 
For the two-mode case we obtain from (14) and (17) for the case of 
low losses 
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L. = 4- (18) 

Finally, we study the steady-state pulse width which follows from (3), 
with i = 1. We also neglect t in this equation, assuming that the pulse 
has spread to a size much longer than the input pulse. From (6), (8), 
(9), (10), and (11) we find 

We used z = L, with L designating the length of the waveguide. The 
factor AT, the width of the pulse in the absence of coupling, is defined as 

at = ik - B L - (2o) 

Compared to the width AT of the uncoupled modes, the pulse length 
in case of coupled modes is improving with length. The pulse width 
formula can again be considered in the two limiting cases. If a 2 — a t <C h 
we have 

at 
At = V§ pp. (21) 

This formula shows clearly that the pulse length shortens with increased 
coupling strength. In the other extreme, h <JC | a 2 — «i | , we have 

At, = 4AT T( ^3n^ (22) 

L(a 2 — a 1 ) Ly 

It appears strange at first that the pulse length now increases as the 
coupling strength is increased. However, in this mode of operation, the 
pulse length is primarily determined by the differential loss of the two 
modes. If both modes travel uncoupled, one will die out while the other 
carries the pulse all by itself. In this case, the pulse width is determined 
only by the dispersion of the surviving mode, which is not included in 
our theory. For h = 0, we thus obtain a vanishing pulse length. As 
the coupling is increased, power is flowing from the lower loss mode to 
the high loss mode so that the pulse width is increased by the different 
delay time of each mode. It is thus clear that a small amount of coupling 
causes the pulse to lengthen. 

Finally, we combine the formula (21) for the low loss case with the 
formula (18) defining the characteristic length that indicates where 
the steady state is reached. We thus obtain the interesting result 
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^--^AT^Y- (23) 

The factor in front of this equation becomes unity if we use k = 4. We 
have seen that this value is large enough to ensure that steady state is 
essentially reached. Equation (23) has been derived by Personick. 4 We 
see from our derivation that the characteristic length l c in Personick's 
formula can be interpreted as the length that is required to reach the 
steady-state power distribution. 

Equation (23) was derived for the case of only two modes. It is 
tempting to use this equation also for the multimode case. In order to 
test the mode dependence of this formula, I solved the four-mode prob- 
lem under the assumption that all off-diagonal elements of h tli vanish 
with the exception of the elements directly adjacent to the main diagonal. 
All non-vanishing elements of h fll were set equal to the same value h. 
For the lossless case, and assuming that the inverses of the group 
velocities are evenly spaced, the following formula was obtained for 
the four-mode case. 

Afc-o-ro-^Arfp 1 )*- (24) 

L. is again defined as the length required to achieve the steady state. 
AT is the length of the uncoupled signal for the four-mode case. Since 
(23) and (24) are essentially identical [(24) is even slightly more favor- 
able], it might be assumed that (23) may hold independently of mode 
number. 

Equation (23) was derived for the case where the coupling coefficient h 
is larger than the loss coefficients a v or a 2 . In the opposite case, where 
the losses determine the rate at which the power distribution approaches 
the steady state, no simple relationship exists between L. and At. 
Formula (23) is quite useful for estimating the length of the Gaussian 
pulse if the distance L. , at which steady state is reached, can be ob- 
served. If it is known that radiation losses are small compared to the 
coupling coefficient h, the conditions exist for which (23) was derived. 
However, it may well be that the region of dominance of radiation 
losses over coupling strength is different for different modes. Experi- 
ments have shown that the coupling mechanism in optical fibers consists 
of two parts. 6 A Rayleigh-type background with a wide mechanical 
Fourier spectrum is responsible for most of the radiation losses, while 
a very sharp peak at zero frequencies of the mechanical power spectrum 
is responsible for most of the coupling between guided modes. The 
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broad spectrum causes more radiation loss for higher-order modes, 
provided that the coupling is caused by core-cladding interface irregu- 
larities. The narrow peak at zero mechanical frequencies couples lower- 
order modes much more strongly than high-order modes, because of 
the closer spacing (in /3-space) of the low-order modes. 3 The combined 
effect of these two spectral regions causes a steady-state distribution 
that favors the lower-order modes. In this situation, it may still be 
possible to estimate the pulse performance of the multimode waveguide 
by ignoring those modes that do not carry power in the steady-state 
distribution and interpret AT as the pulse length that would be obtained 
by the remaining modes in the absence of coupling. The steady-state 
distance L, is best observed by launching only the lowest-order modes 
and measuring the distance that is required until the power versus mode 
number distribution ceases to change its shape. 
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